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Abstract. A new Combinatorial Ricci curvature and Laplacian oper- 
ators for grayscale images are introduced and tested on 2D synthetic, 
natural and medical images. Analogue formulae for voxels are also ob- 
tained. These notions are based upon more general concepts developed 
by R. Forman. Further applications, in particular a fitting Ricci flow, are 
discussed. 



1 Introduction 

Curvature analysis plays a major role in Image Processing, Computer Graphics, 
Computer Vision and their related fields, for many applications, such as recon- 
struction, segmentation and recognition, to list only a few (see, e.g. [3], [10], [16], 
[TT]). Traditionally, the curvature estimation is that of a polygonal (polyhedral) 
mesh, approximating the ideally smooth (C^) image under study, such that the 
curvature measures of the mesh converge to the classical, differential, curvature 
measure of the investigated surface. For surfaces, by far the most important 
curvature is the intrinsic Gaussian (or total) curvature. 

Recently, partly as an offshoot of the great interest generated by G. Perel- 
man's important contribution on the Ricci flow and its application in the proof 
of Thurston's Geometrization Conjecture, and, implicitly of the Poincare Con- 
jecture (see, e.g. ^ for a comprehensive exposition), a flourishing of the study 
of various discrete versions of the Ricci flow (and similar related flows) occurred 
(see |3], 0, m)- 

Ricci curvature measures the defect of the manifold from being locally Eu- 
clidean in various tangential directions. More precisely, it appears in the second 
term of the formula for the (n — l)-volume Q{e) generated within a solid angle 
(i.e. it controls the growth of measured angles) - see Fig. 1. Moreover, 

n - 1 f 

V • Ricci{v) = / K{< V, w >) , 



vol (S^ 



where < v, w > denote the plane spanned by v and w, i.e. Ricci curvature 
represents an average of sectional curvatures. The analogy with mean curvature 
is further emphasized by the following remark: Ricci curvature behaves as the 
Laplacian of the metric g (see, e.g. [2]). It is also important to note that in 
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Fig. 1. Ricci curvature as defect of the manifold from being locally Euclidean 
in various tangential directions (after Berger [l]). Here duo denotes the n- 
dimensional solid angle in the direction of the vector v, and i7(e) the (n — 1)- 
volume generated by geodesies of length e in duj. 



dimension n = 2, that is in the case that is the most relevant for classical Image 
Processing and its related fields, Ricci curvature reduces to sectional (and scalar!) 
curvature, i.e. to the classical Gauss curvature. 

However, both in the more classical context, as well as in the new directions 
mentioned above, smooth surfaces and/or their polygonal approximations con- 
sidered. Unfortunately, smooth surfaces are at best a crude model (and usually 
nothing but a polite fiction), as far as digital and grayscale images, i.e. the stan- 
dard objects of study in Image Processing, are considered, in particular for those 
objects that are not "natural images" , such as images produced using ultrasound 
imaging, MRI or CT. It would be of practical interest to define a proper notion 
of curvature for digital objects, in the spirit of [7J. By "proper" we mean discrete, 
intrinsic to the nature of the spaces under investigation, and not an approxima- 
tion or rough discretization of a differential notion. Moreover, we hope that, by 
doing this, we shall be able to help bridge, at least partially, the divide between 
"Digital" and "Smooth" Image Processing. In addition, since, as we shall show, 
each dimension displays its own Laplacian, we believe our method can produce 
more types of heat fiow, along edges, pixels, etc, thus allowing for better tools for 
the intelligence of higher dimensional data (such as RGB images, images with 
texture, etc.). 

We are fortunate in our quest to be able to rely on the work of R. Forman [5] 
on Combinatorial Ricci curvature and the so called "Bochner Method", where 
he addressed this very problem in the far more general setting of weighted cell 
complexes^ which represent an abstractization of both polygonal meshes and 
weighted graphs. While we succinctly present some of the more general facts 
residing in Forman's work, in this paper we concentrate solely on the case of 
grayscale images with their very special combinatorics and weights, and largely 
defer the study of higher dimensional (color) images and their curvatures and 
Laplacians for further study [15j . 
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The present paper represents an extension of our previous work [14J. More 
precisely, we have augmented the previous version by adding: (a) The formula 
for the Ricci curvature of 3-dimensional cubical complexes (that is, in the case 
of images, for voxels); (b) A discussion on the purely combinatorial version of 
the Ricci curvature and Laplacians; (c) A theoretical comparison of the classical 
and combinatorial Laplacians; (d) A combinatorial diffusion, corresponding to 
upsampling and downsampling (this aspect being augmented by a computational 
example); (e) A discussion of a possible Ricci flow for our combinatorial setting, 
with a concrete suggested direction of study. In addition we have illustrated 
the meaning of Ricci curvature and included more detailed exposition of the 
Algebraic Topology background and ideas. And, last, but not least, we bring 
more (and new) experimental results: on synthetic images (that were not given 
previously), on standard test images and also on medical images. 

2 Forman's Combinatorial Ricci Curvature 

We sketch below some of Forman's main ideas [5|. This requires some technical 
definitions and notations. We do not introduce here the basic technical notions 
in Algebraic Topology and Differential Geometry, and refer the reader to [13] for 
the former and to [2j for the later. 

To generalize the notion of Ricci curvature, in a manner that would include 
weighted cell complexes, one starts from the following form of the Bochner- 
Weitzenbock formula (see, e.g. ^2j) for the Riemann- Laplace operator Hp on p- 
forms on (compact) Riemannian manifolds: 

Bp = dd* + d*d = V;Vp + Curv{R) , (1) 

where V*Vp is the Bochner (or rough) Laplacian and where Curv{R) is an ex- 
pression with linear coefficients of the curvature tensor (Here Vp is the covariant 
derivative operator.) Of course, for cell-complexes one cannot expect such dif- 
ferent iable operators. However, a formal differential exists: in our combinatorial 
context (the operator) "d" being replaced by "9" - the boundary operator of 
the cellular chain complex (see [13j), 

— > Cn{M, R) Cn-i (M, R)-^ ^ Co(M, R) — > , 

were cells are playing in this setting the role of the forms in the classical (i.e. Rie- 
mannian) one. The following definition of the combinatorial Laplacian becomes 
now natural: 

Dp = dd* + d*d : Cp{M, R) Cp{M, R) , (2) 

where 9* : Cp(M, R) Cp+i(M, R) is the adjoint (or coboundary) operator of 
a, defined by: < 9p+iCp+i,Cp > =< Cp+i,5*Cp >p+i , where < •, • >=< •, • >p 
is a (positive definite) inner product on Cp(M, R), i.e. satisfying: (i) < /3 >= 
0, 7^ /3 and (ii) < a, a >= Wa > - the weight of cell a. 
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Forman [5j shows that an analogue of the Bochner-Weitzenbock formula holds 
in this setting, i.e. that there exists a canonical decomposition of the form: 

\Jp = Bp ^ Fp , (3) 

where Bp is a non-negative operator and Fp is a certain diagonal matrix. Bp 
and Fp are called, in analogy with the classical Bochner-Weitzenbock formula, 
the combinatorial Laplacian and combinatorial curvature function, respectively. 
Moreover, if a = is a p-dimensional cell (or p-cell, for short), then we can 
define the curvature function: 

J^p=< Fp{a),a>, (4) 

Fp : Cp ^ Cp being regarded as a linear function on p-chains. For dimension 
p = 1 we obtain, by analogy with classical case, the following definition of discrete 
(weighted) Ricci curvature: 

Definition 1. Let a = be a 1-cell (i.e. an edge). Then the Ricci curvature 
of a is defined as: 

Ric(a) = j^i(a). (5) 

While general weights are possible, making the combinatorial Ricci curvature 
extremely versatile, it turns out (see [5j), that it is possible to restrict oneself 
only to so called standard weights: 

Definition 2. The set of weights {wa} is called a standard set of weights iff 
there exist W\,W2 > such that given a p-cell , the following holds: 

w{a^) = wi • W2 ' 

(Note that the combinatorial weights Wa = ^ represent a set of standard weights, 
with wi = W2 = 1.) Using standard weights, we obtain the following formula for 
polyhedral (and in fact much more general) complexes: 



■ E 



(3P+^>al,(3P+^>aP jP-^ Ka^ ,-fP-^ <aP 



where a < P means that a is a face of and the notation ai \\ a2 signifies that 
the simplices ai and 0^2 are parallel, the notion of parallelism being defined as 
follows: 

Definition 3. Let ai = and a2 = ^'^^ p-cells. ai and a2 are said to 

be parallel (ai \\ a2) iff either: (i) there exists /3 = (5^^'^ , such that ai,a2 < /3; 
or (a) there exists 7 = /3p~^ , such that ol\,ol2 > 7 holds, but not both simulta- 
neously. (For example, in Fig. 1, 61,62,63,64 are all the edges parallel to cq.) 
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Together with the formula above, the (dual) formula for the combinatorial 
Laplacian (see [5 ) is also obtained to be: 



/ ^ -ai,a2,7 / — — 



where eQ,^^Q,2,/35 ^ai,Q!2,7 ^ {"l^+l} depend on the relative orientations of the 
cells. 



3 Combinatorial Ricci Curvature of Images 

Before developing the relevant formulae in the special combinatorial setting of 
the tilling by squares of the plane, as it is usually considered in (Discrete) Image 
Processing, let us first indicate that it is advantageous to use standard weights. 
Such natural weights are proportional to the geometric content (s.a. length and 
area). It follows that the weight of any vertex is w{v) = 0. Bearing this in mind, 
and considering the combinatorics of the square tilling (see Fig. 1), the specific 
form of Combinatorial Ricci curvature for 2D images is: 



Ric(eo) w{eo) 



w{ci) W{C2) J \ W{ci) W{C2) 



• (8) 



For the Laplacian there exists more than one possible choice, depending on 
the dimension p. The simplest, and operating on cells of the same dimensionality 
as the Discrete Ricci curvature, is Di. Because vertices have weight and adja- 
cent cells have opposite orientations. Equation (7) becomes, in this case (using 
the notation of Fig. 1): 

w{ci) W[C2) 

The formula for the Combinatorial Bochner Laplacian follows immediately: 

5i(eo) = 01(60) -Ric(eo). (10) 

Instead of computing a Laplacian along the edge eo, one can compute a 
Laplacian operating across the edge, namely 02(61,02). Since no 3-dimensional 
cells exist, the first sum in Equation (7) vanishes. Hence, we have (up to sign): 



(11) 
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Fig. 2. The elements appearing in the computation of the Combinatorial Ricci 
curvature of edge eo in the geometry of pixels. 



Remark 1. By "weighing" the combinatorial formula for the curvature function 
(cf, [5j) we obtain: 



(12) 



(Here denotes the number of elements of the set X.) For p = dim a = 1, a 
simplified, purely combinatorial version of Formula (|8]) is also obtained: 

Ric(eo) = w{ci) + w{c2) — w{ei) - w{e2) - w{e3) — w{e4) + 2. (13) 

However, the above combinatorial version of the Ricci curvature does not provide 
considerable geometric insight and, therefore, does not yield interesting results. 

In the case of the cubical geometric configuration present in Digital Images, 
Formula [6] becomes: 



Ric(eo) w{eo) 



E 



^(eo) 



E 



<c'>eo ^ ^ c0<e2 ^ ^ 



E 

e||eo,e/eo 



E 



c^>e,c^>eo 



^"^(60)^(6) 

w{c'^) 



E 

cO<e,cO<eo 



^(cO) 



Since, as we have already noted, for digital images the vertices' weights are 
always 0, we obtain the following expression for Ric(eo): 



Ric(eo) ^(eo) 



1 



7^ ^(Ci)^ 



w{ci) 



(14) 



(see Fig. 2, right). 
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4 Experimental Results 

Before commencing any experiments with the combinatorial Ricci curvature in 
the context of images, we had to choose a set of weights for the 2-, 1- and 0- 
dimensional cells of an image, that is for squares (pixels), their common edges 
and the vertices of the tilling of the image by the pixels. Any such choice should, 
obviously, be as natural and expressive as possible for image analysis. The choice 
of weights was motivated by two factors: the context of Image Processing, where 
a natural choice for weights imposes itself (see below) and the desire (and, indeed, 
sufficiency, see Section 2) to employ solely standard weights. 

Since natural weights have to be proportional to the dimension of the cell, it 
follows immediately that the weight of any vertex (0-cell) has to be 0. Moreover, 
in the beginning, it is natural to choose wi = 1 and W2 = length of a cell. A 
somewhat less arbitrary choice for the length (i.e. basic weight) of an edge, 
would be Length(e) = (dimension of the picture)"^, hence that for the area (i.e. 
basic weight) of a pixel a being Area(a) = (dimension of the picture)"^. The 
proper weight for a cell a should, however, take into account the gray-scale level 
(or height) of the pixel in question, i.e. = ■ Area(a). This will become, 
so we hope, clearer in the following paragraph. The natural weight for an edge e 
common to the pixels a and (3 is \ha — hp\. {A less "purely" combinatorial choice 
of cells and weights is discussed in [15j.) 

Note that, given an edge e, the Ricci curvature Ric(e) represents in a way a 
generalized mean of the weights the cells parallel to e. Therefore, it represents 
a measure of flow in the direction transversal to e. It follows, that, contrary 
perhaps to intuition, this type of Ricci curvature (and the Bochner Laplacian 
associated to it) in direction, say, parallel to the x-axis, is suitable for the detec- 
tion of edges and ridges in the ^-direction. On the other hand, since scalar (i.e. 
Gauss) curvature, is associated to each pixel, that is to each square of the tes- 
sellation, to compute the Gaussian curvature one has to compute the arithmetic 
mean of the Ricci curvatures of edges of the square under consideration - see 
Fig. 2. (A similar argument holds if one wishes to compute the 1-Laplacians, Di 
and 5i, of a given pixel.) The difference between the Ricci curvature computed 
in the horizontal and vertical directions, as well the "true", i.e. average Ricci 
curvature can be seen in Fig. 3 and Fig. 4. As Fig. 5 illustrates, the Combina- 
torial Ricci curvature we introduced herein allows even for a non-optimal choice 
of weights, a very good approximation of Gaussian curvature of surfaces (i.e. 
for gray-scale images). Here, classical Gaussian curvature was computed using 
finite element methods standard in Image Processing - see [16 . In contrast, 
both the Bochner (and Riemann) Laplacian sharply diverge from the classical 
one, e.g. the one obtained by using the standard Matlab function (see Fig. 6). 
This is not too surprising, given the fact that such a comparison is, in a sense, 
not relevant, due to the different dimensionality of the two concepts: The Com- 
binatorial and Bochner Laplacian are, as stressed above, associated to edges, 
hence 1-dimensional (this being underlined by the notation: □i(eo) and 5i(eo), 
respectively). In contrast, the classical Laplacian is a pointwise function, (and, 
in its discrete setting, associated to the vertices of a mesh) hence 0-dimensional. 
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Fig. 3. Computing the Ricci curvature for a pixel: Synthetic test image (left), 
Ricci curvature computed in horizontal (middle) and vertical (right) directions. 



However, the Bochner Laplacian proves to be an excellent detector of "sharp" 
edges (see Fig. 7), therefore it may prove to be useful for contour detection and 
for segmentation. A "combinatorial diffusion" was also performed, either by 
upsampling, that is by division of the squares (pixels) into sub-squares of length 
1/2 or 1/3, or by downsampling, i.e. by the fusion of 4 or 9 squares into larger 
ones. However, using the fusion of 9 squares produces more inferior results.) The 
weight after such upsampling would be, for instance ha/^- The effect of this last 
process on 01(60) and on 5i(eo) can be seen in Fig. 8. 



5 Future Work 

We briefly discuss below some of the natural and/or seemingly required direc- 
tions of further study: 

1. Evidently, the first task in testing the efficiency of the combinatorial Ricci 
curvature and Laplacian in medical imaging is to experiment with voxels, 
that is to apply the apparatus introduced herein to the analysis of volumetric 
data. Such 3- (and even 4-) dimensional manifolds and their evolution in time 
is most relevant, for instance, in the analysis of cardiac MRI. This brings us 
to the following point: 

2. As already mentioned in the introduction, the full power of the Ricci cur- 
vature is revealed in the general heat-type diffusion and in discrete versions 
of the Ricci flow (and other related flows), that were introduced and ex- 
perimented with elsewhere (0, [6], [8], [H])). It is only natural to strive to 
develop and experiment with a discrete version of the Ricci flow correspond- 
ing to the combinatorial Ricci curvature introduced herein. The Ricci flow 
is given by 

|| = -2Ric(i;(i)), 

where g denotes the metric of the manifold, (see e.g. [I2]). Various discretiza- 
tions are possible for this flow. However, it seems that given the fact that 
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Fig. 4. Computing the Ricci curvature for a pixel: The "Camera Man" test 
image (top, left), Ricci curvature computed in horizontal (top, right) and vertical 
(bottom, left) directions, and the Ricci (averaged) curvature (bottom, right). 



the Combinatorial Ricci curvature is an edge measure, the most natural 
discretization in our context is the following one, inspired by [llj (see also 

my- 

^ = -2Ric(e,,), 

where lij denotes the length of the edge Cij . Experimental work on this type 
of flow is in progress [15j. 
3. As already noted, while we prefer, both for theoretical as well as for practical 
reasons, to work with standard weights, Forman's combinatorial version of 
Ricci curvature is extremely versatile. Even if restricting oneself to using 
standard weights, i.e. proportional to the p-dimensional geometric content 
(p- volume), one still has freedom in choosing the weights W2 and especially 
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Fig. 5. Gauss (left) and Combinatorial Ricci (right) comparison. Note the edge 
detecting capability if the Combinatorial Ricci curvature. 



wi (see Section 2). Hence, to obtain best result, one can experiment in order 
to empirically determine, by using, e.g. variational methods, the optimal 
standard weights for a given application of the method. 
4. The Combinatorial Laplacian is closely connected (by its very definition) 
to the cohomology groups of the cellular complex on which it operates (see 
[5]). It is natural, therefore, to apply the results and methods of [5] for the 
estimation of the dimension, and in some cases even the computation, of the 
cohomology groups (and by duality, of homology groups) of images. (See |9] 
on this direction of study in Image Processing.) 
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Fig. 7. Axial brain scan image (top, left), its Combinatorial Ricci curvature 
(top, right), the Combinatorial (Forman) Laplacian □i(eo) (bottom, left), and 
the Bochner (rough) Laplacian Bi{eo) (bottom, right). (Recall that Bi{eo) = 
□i(eo)-Ric(eo).) 
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Fig. 8. Comparison of Different Laplacians after upsampling, by fusion of 4 
adjacent squares: From left to right: The Combinatorial (For man) Laplacian 
□i(eo) (left), and the Bochner (rough) Laplacian Bi{eo) (right). 



